Exchange-correlation kernels for excited states in solids 
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The performance of several common approximations for the exchange-correlation kernel within 
time-dependent density-functional theory is tested for elementary excitations in the homogeneous 
electron gas. Although the adiabatic local-density approximation gives a reasonably good account of 
the plasmon dispersion, systematic errors are pointed out and traced to the neglect of the wavevector 
dependence. Kernels optimized for atoms are found to perform poorly in extended systems due to an 
incorrect behavior in the long- wavelength limit, leading to quantitative deviations that significantly 
exceed the experimental error bars for the plasmon dispersion in the alkali metals. Copyright 2000 
by The American Physical Society. 



PACS numbers: 71.45. Gm, 71.15.Mb, 71.10.Ca 
I. INTRODUCTION 

Exchange and correlation effects are crucial for under- 
standing the properties of interacting many-electron sys- 
tems but notoriously difficult to implement accurately 
in ah initio computational schemes. In^the Kohn-Sham 
formulation of density-functional theoryEJ they are incor- 
porated into the exchange-correlation potential, which is 
a functional of the electron density but in practice is not 
known exactly except for simple model systems. Hence 
approximations such as the local-density or generalized 
gradient approximations are needed. In many cases these 
yield accurate total energies and related ground-state 
quantitiesB However, the development of experimental 
devices that allow, in principle, to track a single electron 
and the emergence of new fields such as surface pho- 
tochemistry place increasing emphasis on the study of 
excited states. Unfortunately its variational foundation 
prevents a straightforward application of the Kohn-Sham 
scheme to electronic excitations, for which numerically 
expensive Green function techniques or, in case of small 
systems, quantum-chemical methods were traditionally 
employed. p. 

Time-dependent density-functional theoryQ promises 
an appealing alternative. Originally designed to explore 
time-dependent phenomena, it was recently realized that 
it can also be exploited to investigate optical excitations, 
which involve the creation of electron-hole pairs. □ These 
are to be distinguished from photocmission states, which 
can be determined by calculating the self-energy correc- 
tion to the quasiparticle [band structure, for instance in 
the GW approximation.d The optical excitation ener- 
gies correspond to the poles of the full linear density- 
response function, which can, in principle, be obtained 
exactly within time-dependent density-functional theory. 
This approach has since been applied to excitjed states 
not only in the context of quantum cheniifiteylj'Bl3 but 
also, more recently, in solid-state physics.ciy The pro- 
cedure starts from time-independent ground-state prop- 
erties, such as the Kohn-Sham orbitals and eigenvalues, 
which are conveniently obtained by conventional means. 
However, besides static exchange and correlation embod- 



ied in the exchange-correlation potential, dynamic effects 
due to the time-dependent perturbation must also be ac- 
counted for. The latter are described by the so-called 
exchange-correlation kernel. Formally the kernel is a 
functional derivative of the exchange-correlation poten- 
tial, evaluated at the unperturbed ground-state density, 
but as the exact potential is unknown and must in prac- 
tice be approximated by a parameterization, independent 
expressions are often used for both quantities. 

Chemical studies have repeatedly emphasized that an 
accurate description of the exchange-correlation poten- 
tial is particularly important Joij— excited-state calcula- 
tions of atoms and molecules.EJ^E^ Although the local- 
density and generalized gradient approximations often 
yield good total energies, the corresponding potentials 
fail to capture the correct Coulomb-like asymptotic be- 
havior and instead decay exponentially. As a result, 
many unoccupied states that really should be bound are 
pushed to higher energies and merge with the contin- 
uum. The accessible excitation spectrum is then poorly 
rendered. Much effort is therefore invested into better ex- 
pressions for the exchange-correlation potential. In com- 
parison, the kernel is considered to be less important, and 
the adiabatic local-density approximation is often chosen 
for convenience. 

It is not clear whether these findings apply equally to 
solids. First, the problem of possible unbound states 
becomes irrelevant in bulk materials. Second, nonlocal 
dynamic exchange and correlation effects, which are ne- 
glected in the adiabatic local-density approximation, nat- 
urally become more prominent in extended systems with 
delocalized wavef unctions. For solids the errors intro- 
duced by approximations for the potential and the ker- 
nel could therefore be of similar magnitude. However, al- 
though several new parameterizaticm^Qr. the kernel were 
recently proposed in the literatureJjll3'EZl these have not 
yet been systematically applied to solids, despite first 
encouraging attempts B In fact, almost all calculations 
reported so far, which fotps either on the dielectric re- 
sponsQ O k f s emiconductarsOilj or the plasmon dispersion in 
simplcOEj and nobleEj metals, employ the local-density 
approximation both in the potential and the kernel. In 
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contrast to atoms, this approach seems to work reason- 
ably well for solids, but improvements are still desirable. 

In order to understand and quantify the error intro- 
duced by approximations to the kernel, more detailed, 
studies for solids are necessary. Lein and co-workeralS 
have recently compared the correlation energy of the ho- 
mogeneous electron gas for different parameterizations, 
but the energy is an integrated quantity that principally 
reflects average weight distributions and is less sensitive 
to variations in the small-scale structure of the kernel. 
However, the latter have a significant influence on the 
excitation spectrum, which is given by the position of 
poles in the linear density-response function. In this pa- 
per we therefore concentrate on the dispersion relations 
for the plasmon frequency and lifetime, assessing the per- 
formance of several kernels currently circulated in the lit- 
erature. By applying them to the homogeneous electron 
gas, where the exchange-correlation potential is a trivial 
constant, we are able to isolate the error due to the kernel 
and make a systematic comparison. Indeed, we find that 
the choice of parameterization plays an important role, 
and that inappropriate kernels optimized for atoms give 
rise to quantitative deviations that significantly exceed 
the experimental error bars for the plasmon dispersion 
in the alkali metals. We believe that the findings pre- 
sented in this paper can not only be extended to plasmons 
in real materials but generally apply to excited states 
in solids that, like plasmons, are based on charge rear- 
rangements. A prominent example are charge-transfer 
excitations in surface-adsorjaate systems that occur dur- 
ing photoinduced reactionsJlj 

This paper is organized as follows. In Sec. || we give 
an outline of our computational method. The kernels 

and in Sec. |^ we 
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considered here are listed in Sec. 
present the numerical results together with a discussion. 
Finally, in Sec. we summarize our conclusions. Ryd- 
berg atomic units are used throughout. 



II. COMPUTATIONAL METHOD 

Within linear response theory the true many-body 
density-density response function is defined as 



X(r,r';t-i') 



Sn{r, t) 



(1) 



where Sn{r,t) indicates the density change induced by 
an external perturbation (5F°^*(r', t'), and the functional 
derivative is evaluated at the static external potential 
corresponding to the unperturbed ground-state density. 
Likewise, the Kohn-Sham susceptibility 



X°(r,r';i-t') = 



6n{r, t) 
(5F°ft(r',i') ' 



(2) 



describes the response of the associated noninteracting 
system with the same electron density due to a change 



SV^iy'^t') = 5V'"^*(r',t') + 5V"{T',t') + SV^ir'^t') in 
the effective potential, which includes the Hartree and 
exchange-correlation contributions. It can be calculated 
explicitly in frequency space according to 



X"(r,r';a;) = 2^(/.,-/.) 



^.(r)^:,(r)^*(r')^.'(r') 



UJ 



(3) 



from the static Kohn-Sham orbitals Lp^, and eigenvalues 
Cij. The symbol fn indicates the Fermi occupation num- 
bers and 77 is a positive infinitesimal. The true density- 
density response function is obtained by relating it to 
the Kohn-Sham susceptibility through the chain rule for 
functional derivatives, which may be written in the form 



d\" Miy, r"; c.)x(r", r'; - x°(r, r'; u:) 



(4) 



The operator 



M(r, r'; lo) = S{r -^')- J d^^" X°(r, r"; 

is related to the dielectric function and the exchange- 
correlation kernel is defined as 



r{r,v';t-t') = 



(6) 



where the functional derivative is evaluated at the un- 
perturbed ground-state density. 

The many-body density-density response function has 
poles at the exact excitation energies of the interacting 
electron system. On the other hand, x" has poles at 
the excitation energies of the corresponding Kohn-Sham 
system, which are in general different. Hence the singu- 
larities of X in Eq. must be cancelled by zeroes of the 
operator M . Independent of the system characteristics 
and the nature of the excited states, this provides a con- 
venient starting point for calculating the excitation spec- 
trum. For atoms it is possible to expand all quantities in 
a Laurent series around a particular Kohn-Sham energy 
difference en—eni , which leads to an explicit expression for 
low-order corrections to the transition energy.a This pro- 
cedure is not appropriate for solids, however, because in 
infinite systems the Kohn-Sham transition energies form 
a continuumEj Although the resulting structure of 
can still be described in terms of effective poles, these 
are located off the real axis and no longer correspond to 
individual energy differences in the denominator of Eq. 
(^. For the homogeneous electron gas we therefore de- 
termine the plasmon dispersion 0.{q) by a direct search 
for the zeroes of 



M(g,^) = 1 - x"(9,^) Hq) + riq^oj)] 



(7) 
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FIG. 1. At the border of the electron-hole pair continuum, 
indicated by the shaded region, a crossing between the plas- 
mon dispersion and a satellite excitation occurs. 



in reciprocal space, which is also the most accurate pro- 
cedure. The Fourier transform of the Coulomb potential 
is v{q) = 47r/(j'^, and the Kohn-Sham susceptibility x'Jjs 
given analytically by the dynamic Lindhard function.Lj 
In our implementation we generally calculate the ze- 
roes of M in the complex frequency plane, which allows 
us to obtain both the plasmon dispersion and the corre- 
sponding lifetimes. This approach requires the simulta- 
neous solution of a set of two nonlinear equations for the 
real and the imaginary part, for which we employ an iter- 
ative procedure. At g = the classical plasma frequency 
LOp = {iTTny^"^ forms a convenient starting point, while 
for finite q the previously calculated solution for a smaller 
wavevector is utilized. We have confirmed that the itera- 
tion is stable and convergent. Special attention must be 
paid at ^}{q) ~ ^q^ + qkp, however, where the structure 
of the excitation spectrum becomes very complicated and 
a bifurcation occurs. This bifurcation, illustrated in Fig. 
|l| for the random-phase approximation and — 4, where 
rs denotes the Wigner-Seitz radius, is due to the cross- 
ing of the plasmon dispersion with a satellite excitation. 
In order to discriminate between the principal resonance 
and its satellite, we examine the spectral function 



S{q,uj) 
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nv{q) M{q,u}) 



(8) 



on the real frequency axis, where the plasmon peak can 
easily be identified as the dominant feature. 

The imaginary part ofn{q) reflects the finite lifetime of 
plasmons in the homogeneous electron gas. The most im- 
portant decay mechanism is scattering into electron-hole 
pairs, which dominates whenever energy and momentum 
conservation allow the promotion of an electronpLnto a 
previously unoccupied state above the Fermi level.o This 
region, bounded by the two lines -^q^ ± qkp, where kp 
denotes the Fermi wavevector, is shown shaded in Fig. ||. 



Real solutions for the plasmon energy, suggesting an infi- 
nite lifetime, only exist for frequency-independent kernels 
outside the electron-hole pair continuum. This behavior, 
well known from the random-phase approximation, con- 
stitutes a physically implausible artefact that stems from 
the neglect of more intricate decay channels. 



III. EXCHANGE-CORRELATION KERNELS 

In this paper we have considered the following approx- 
imations for the exchange-correlation kernel. 

(a) In the random-phase approximation (RPA) all dy- 
namic exchange-correlation effects are ignored by setting 



RPA 







(9) 



(b) The adiabatic local-density approximationcS 
(ALDA) equals the long-wavelength limit 



ALDA 



0) 



(10) 



of the static exchange-correlation kernel of the homoge- 
neous electron gas. It is readily expressed in terms of the 
exchange-correlation energy per particle Chom 



ALDA 



(11) 



Owing to its computational simplicity, the ALDA has 
become the standard approximation in time-dependent 
density-functional theory. It has already been employed 
in calculations of the plasmon dispersion for solids.EJ~Ej 
(c) In their original application of time-dependent 
density-functional theary to excited states, Petersilka, 
Gossmann, and GrossQ (PGG) derived an approximate 
exchange-only kernel in the spirit of the optimized ef- 
fective potential method.ES This approach has the ad- 
vantage that the corresponding exchange potential has 
the proper Goulomb decay. The kernel is constructed 
from the Kohn-Sham orbitals and hence only depends 
implicitly on the density. Designed for small atoms, the 
PGG formula is identical to the exact exchange kernel 
for two-electron systems, but deviations are expected for 
extended systems. In particular, the frequency depen- 
dence of the exact exchange kernel, which, in principle, 
can also be calculatedJSil is ignored. In momentum space 
the PGG kernel is given hyt3 



/pgg(9) 



3n 

Tofcf 



ii + 2g2 



(2Q4 - 10Q2) hi 



lOQ In 



1 + Q 

11-01 

(12) 



with Q = q/2kp. _ 

(d) Burke, Petersilka, and GrossEB (BPG) recently pro- 
posed a hybrid formula that was shown to improve the ex- 
citation spectra of small atoms. It combines expressions 
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for symmetric and antisymmetric spin orientations from 
different approximations in a spin density-functional for- 
malism. For the unpolarized homogeneous electron gas 
this kernel reduces to 



/bpg(9) 



[/pgg,tt('?) 



ALDA.TiJ 



(13) 



(e) An essentially exact parameterization of the static 
exchange-correlation kernel for the homogeneous elec- 
tron gas was given by Corradini, Del Sole, Onida, and 
PalummoO (CDOP), who used thx^Montc Carlo results 
of Moroni, Ceperley, and SenatoreE3 for the static local- 
field factor G{q) = —f^'^{q)/v{q). Unlike the original 
data, the fit 



/cDOp(9) 



An 



CQ^ 



BQ^ 



aQ e 



(14) 



with Q = q/kp is not restricted to metallic densities, 
because it incorporates the known asymptotic limits for 
high and low densities. The parameters a, /3, S, C, and 
g depend on Vg and are listed in Ref. |l^. By construction, 
the CDOP kernel becomes identical to the ALDA in the 
long- wavelength limit. 

(f) Finally we consider a parameterization of the dy- 
namic local-field factor of the homogeneous electron gas 
proposed by Richardson and AshcroftlHj (RA), includ- 
ing the corrections given in Ref. |l^, which stems from 
the summation of self-energy, exchange, and fluctuation 
terms in the diagrammatic expansion of the polarization 
function. It satisfies many important sum rules and in- 
cludes exact asymptotic expressions for small and large 
wavevectors. At intermediate wavevectors and frequen- 
cies it provides a realistic description of the position and 
magnitude of extrema, which are related to the pair dis- 
tribution function evaluated at zero separation. Because 
of this careful derivation we believe the RA expression to 
be very close to the exact dynamic exchange-correlation 
kernel of the homogeneous electron gas and give an ac- 
curate account of the plasmon dispersion. In the absence 
of experimental data we therefore use the RA results as 
a reference in order to assess the performance of simpler 
approximations. The parameterization 

/rK?, ^) = -% [Gs(Q, U) + G„(Q, [/)] (15) 

r 

with Q = ql2kp and U — uj/Akp was originally given on 
the imaginary frequency axis, so we use a continuation to 
the full complex plane. The local-field factor Gs describes 
screened exchange, fluctuation, and self-energy effects in 
the irreducible polarizability, while Gn accounts for the 
change in occupation numbers due to correlation. The 
parameterized forms of both are listed in Ref. 

Whenever the exchange-correlation energy is needed 
as an input, we use the parameterization by Perdew 
and W^gp3 of the Monte Carlo data by Ceperley and 
Alder .Ea The pair-correlation function that enters the RA 
kernel is taken from the same authors.E3 



IV. RESULTS AND DISCUSSION 

Before presenting our numerical results we first discuss 
what can be deduced from an analytic expansion of the 
plasmon dispersion il{q). By solving Mig,Q{q)) = up 
to second order in q we obtain the seriesB 



fl{q) = LOp 
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10/L^j^p 



8tt 



0{q^) 



(16) 



where fcxp indicates the Thomas-Fermi wavevector. As 
long as /'"^ does not diverge, all curves should approach 
the classical plasma frequency in the long-wavelength 
limit. The kernel only introduces corrections beyond 
the RPA in second order, where the element f^'^{Q,ujp) 
appears. The ALDA evidently contains the right long- 
wavelength limit of the kernel, but its neglect of the 
frequency dependence still introduces an error in the 
parabolic term. By construction, the CDOP formula 
becomes identical to the ALDA in the long-wavelength 
limit and thus produces the same second-order term. Be- 
ing static approximations, the PGG and BPG function- 
al are also evaluated at w = rather than the plasma 
frequency. However, they do not approach the correct 
long-wavelength limit of the homogeneous electron gas 
and therefore generate a different parabolic coefficient. 
The RA kernel, which incorporates the full frequency de- 
pendence, is the only parameterization that is formally 
exact beyond the trivial zeroth order. 

The numerically calculated plasmon dispersions for 
rg — 4: are shown in Fig. ||. As predicted, all curves start 
at the classical plasma frequency. For small wavevectors 
only a small spread of the results is observed, because 
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FIG. 2. Plasmon dispersion for the homogeneous electron 
gas at r-j, = 4, calculated with different approximations for 
the exchange-correlation kernel. The electron-hole pair con- 
tinuum and the resulting nonzero imaginary part of the plas- 
mon frequency in this regime are also marked. 
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the factor 9/lOfc^p in Eq. ( p^ in general outweighs the 
contribution of the kernel. However, a slight downward 
shift compared to the RPA is clearly visible for all non- 
trivial approximations, because dynamic exchange and 
correlation effects combine to lower the energy of the 
electron system. The ALDA and the CDOP formula pro- 
duce curves that are initially very close to the RA result 
we use for reference, indicating that the neglected fre- 
quency dependence is of little consequence as long as the 
correct long-wavelength limit is reproduced. This point 
is emphasized by the relatively large deviation for the 
static PGG kernel, which stems precisely from its incor- 
rect behavior at g ^ 0. The BPG curve, as expected, 
lies between the ALDA and PGG results. 

To demonstrate that these observations are representa- 
tive, in Fig. ^ we show the behavior of the plasmon energy 
over a large density range. The curves are calculated for 
q = 0.5(7c, where qc indicates the critical wavevector cor- 
responding to the onset of damping due to electron-hole 
pair excitations in the RPA. Note that the results are 
scaled in units of the Fermi energy ep^ which is itself a 
function of rg. The RPA and ALDA curves are practi- 
cally indistinguishable on the scale of the figure. In the 
high-density limit all parameterizations tend to the RPA 
result, which is the correct trend. The deviation between 
the RA dispersion and the other curves increases approx- 
imately linearly with r^. 

At larger wavevectors, where the parabolic expansion 
( p^ ) is no longer valid, the differences between the con- 
sidered approximations become more pronounced. The 
dispersion resulting from the static ALDA kernel remains 
close to the RPA at too high energies, while the CDOP 
result begins to deviate slightly from the RA curve after 
the onset of damping in the electron-hole pair continuum. 
This discrepancy must be attributed to the static nature 
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FIG. 3. Behavior of the plasmon energy as a function of 
the density evaluated at g = O.Sqc, where qc indicates the 
critical wavevector corresponding to the onset of damping due 
to electron-hole pair excitations in the RPA. 
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FIG. 4. Behavior of the plasmon energy in the region of 
electron-hole pair excitations at g = 1.2gc- 



of the CDOP kernel. Furthermore, it can be seen that 
the strong downward shift of the exchange-only PGG for- 
mula leads to an even larger error in absolute terms than 
the underestimation of dynamic exchange and correla- 
tion effects in both the ALDA and RPA. The hybrid 
BPG formula, which combines the PGG and ALDA pa- 
rameterizations, profits from a partial cancellation of er- 
rors but improves only marginally upon PGG. In Fig. ^ 
we again show the plasmon energy as a function of the 
density for q = 1.2gc- 

Due to decay into electron-hole pairs in the damped 
regime, the plasmon energy acquires a nonzero imaginary 
part, also displayed in Fig. ^, whose inverse is propor- 
tional to the lifetime. As a general rule we find that all 
kernels yield the same quality of approximation for the 
imaginary part as they do for the real part of the plas- 
mon energy. At small wavevectors, as discussed above, 
static kernels predict a vanishing imaginary part, which 
corresponds to an unphysical infinite lifetime. This arte- 
fact results from modelling f^'^ as a purely real quantity 
by evaluating it at w = 0. In fact, the exact kernel has 
a finite imaginary part at nonzero frequencies, which for 
small wavevectors is related to thajnulti-pair component 
of the susceptibility according toEJ 



Imr=(g,a;) 



[v{q)flu,x^^{q.^) 



(17) 



Such multi-pair decay channels are ignored in the RPA 
and related schemes, which is ultimately the reason for 
their qualitatively wrong behavior. Mermin's modifica- 
tion of the Lindhard dielectric function avoids the prob- 
lem of infinite lifetimesj^ but the correction based on re- 
laxation times is introduced in a phenomenological man- 
ner that makes it unsuitable for ah initio calculations. In 
this study only the dynamic RA parameterization cor- 
rectly predicts a finite plasmon lifetime over the entire 
frequency range. However, outside the electron-hole pair 
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FIG. 5. Calculated plasmon dispersion for = 4.0 com- 
pared to experimental data for sodium from electron en- 
ergy-loss spectroscopy (Ref. ^2|). The theoretical curves have 
been shifted rigidly to the experimental value at g = in 
order to account for core-polarization effects not included in 
this electron-gas treatment. 



wavevectors this parameterization incorrectly predicts a 
negative dispersion that is not observed in Na, although 
negative dispersion does occur in heavier alkali metals, 
such as Cs.Ej The reasons for this anomalous behavior 
are still controversial.EdS Obviously, in such situations 
a poor parameterization may become a serious obstacle 
for theoretical interpretations. After the onset of damp- 
ing due to electron-hole pair excitations, the experimen- 
tal dispersion flattens slightly. As shown in Fig. |^, the 
theoretical results exhibit the same effect, but the un- 
shifted RA and ALDA curves only cross the border of 
the damped regime at larger critical wavevectors. Hence 
in Fig. H quantitative agreement cannot be expected for 
large wavevectors due to the different physical situations. 
The PGG curve, on the other hand, lies below the RA 
result and consequently enters the damped regime at a 
smaller critical wavevector, but the good agreement with 
the experimental dispersion for large wavevectors in Fig. 
H is clearly fortutious. 

V. SUMMARY 



continuum the imaginary part of the plasmon energy is 
several orders of magnitude smaller than the real part 
and hence not discernible in the plot. 

The good agreement between the static CDOP param- 
eterization on the one hand and the dynamic RA result 
on the other over a large wavevector and density inter- 
val indicates that the frequency dependence of the ker- 
nel plays a weak role for the plasmon dispersion. In 
contrast, the significant discrepancy between static ap- 
proximations like ALDA that contain the correct long- 
wavelength limit and others such as PGG, which do not, 
suggests that a correct parameterization of the wavevec- 
tor dependence is crucial. Similar conclusions concerning 
the relative importance of the frequency and wavevector 
dependence_sYere recently also reported for the correla- 
tion energy.Lfl This was not eatirely surprising, however, 
because a frequency analysiala reveals that the dominant 
contribution to the energy in any case comes from the 
low-frequency limit, which is, by design, contained cor- 
rectly in all static approximations. 

To emphasize the error that may arise from an inap- 
propriate kernel, in Fig. |5| we compare theoretical re- 
sults for Ts = 4.0 with experimental data for sodium 
from electron energy-loss spectroscopy.!^ The theoreti- 
cal curves have been shifted rigidly to the experimental 
value at g = in order to account for core-polariza±ion 
effects not included in this electron-gas treatment .lj At 
small wavevectors the RA result is in excellent agreement 
with the experimental dispersion. Although the ALDA 
correctly reproduces the qualitative features, its grow- 
ing deviation from the theoretical reference curve soon 
exceeds the experimental error bars and becomes quite 
pronounced at intermediate wavevectors. This is even 
more obvious in the case of PGG. Furthermore, for small 



In this paper we have tested several common approxi- 
mations for the exchange-correlation kernel by examining 
the plasmon dispersion of the homogeneous electron gas. 
First of all, we have found that the influence of the ker- 
nel is indeed significant, giving rise to large differences 
between the calculated dispersion curves. The ALDA 
performs reasonably well, although it underestimates dy- 
namic exchange-correlation effects embodied in the ker- 
nel and improves only little upon the RPA. A better 
quantitative scheme is therefore desirable. In this respect 
our results, in particular the good agreement between the 
dynamic RA parameterization and the static CDOP ker- 
nel, both of which are presumed to be very accurate, 
suggest that the neglect of the frequency dependence is 
of little consequence, giving rise to small deviations only 
at large wavevectors in the electron-hole pair continuum. 
The error of the ALDA thus stems largely from its local 
nature, and extensions should focus on a better descrip- 
tion of the wavevector dependence. The challenge of this 
task is illustrated by the fact that some of the explicitly 
nonlocal parameterizations we considered, notably PGG, 
actually lead to worse results although they are known to 
improve excitation spectra in small atoms. This appar- 
ent paradox may be understood by the significance of the 
long-wavelength limit for the homogeneous electron gas, 
which determines the leading order of the plasmon dis- 
persion and is contained correctly in the ALDA but not 
in the PGG kernel. In localized systems such as atoms, 
on the other hand, the long-wavelength limit is less rele- 
vant, whereas dynamic exchange effects contained in the 
PGG kernel may play an important role. This lack of 
transferability should encourage specific approximations 
for solids. The CDOP kernel, which derives from the ho- 
mogeneous electron gas, seems a step into the right direc- 
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tion, although its performance for real materials has not 
been fully explored yet. Finally, we demonstrated that 
exchange-correlation kernels optimized for small atoms 
may lead to quantitative, and occasionally qualitative, 
deviations in the plasmon dispersion for solids that sig- 
nificantly exceed the corresponding experimental error 
bars and may affect theoretical interpretations. 
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